Two-dimensional arrays of superconducting strips as dc magnetic metamaterials 
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We theoretically investigate the magnetic response of two-dimensional arrays of superconducting 
strips, which are regarded as essential structures of dc magnetic metamaterials. We analytically 
obtain local distributions of the magnetic field for the ideal complete shielding state (i.e., A/io — >■ 0, 
where 2w is the strip width, A — /d is the Pearl length, A is the London penetration depth, and 
d is the strip thickness), and derive effective permeability by averaging the local field distributions. 
We also perform numerical calculations for a realistic case, taking finite A/w > into account. 
We investigate two types of strip arrays: a rectangular array and a hexagonal array. The resulting 
effective permeability has large anisotropy that depends on the dimensions and arrangement of the 
superconducting strips, and the hexagonal array is found to be more advantageous for obtaining 
large anisotropy than the rectangular array. 

PACS numbers: 74.25. Ha, 74.78. Fk, 81.05.Xj 



I. INTRODUCTION 

The London model describes the macroscopic electro- 
magnetic behavior of a superconductor when fields and 
current are weak and vary on a length scale A (the London 
penetration depth) li In many cases, the superconducting 
sample has large dimensions such that the currents flow 
only along the surface of the sample, which is completely 
shielded from the external magnetic field. This complete 
shielding state corresponds to the limit A — >■ 0. Another 
macroscopic model for describing superconductors with 
strong vortex pinning is the critical state model. ^ This 
model assumes that the current density induced inside a 
superconductor cannot exceed the critical current density 
Jc- The complete shielding state is also observed at the 
high Jc (or low applied field) limit in the critical statei^"— 

In addition to the trivial cases of infinite cylinders or 
slabs in a longitudinal applied field, distributions of sur- 
face current density in flat superconducting samples un- 
der uniform transverse magnetic fields can be analytically 
determined for the complete shielding state in several ge- 
ometries, such as thin disks or strips,^ coplanar pairs of 
thin strips,^ periodically arranged thin stripS)^ and thick 
infinite tapes^. In Ref. |7|, the exact analytical expres- 
sions were presented for the sheet current density and 
magnetic field profiles in an infinite stack and an infinite 
coplanar array of thin superconducting strips in the com- 
plete shielding state. Numerically, the sheet currents can 
also be calculated for thin plates of arbitrary shape in a 
perpendicular applied field. ^"^^ 

Within the London model taking finite A > into ac- 
count, complete shielding is not obtained because the 
magnetic field penetrates the superconductor. In a su- 
perconducting strip of thickness d and width 2w (where 
d <^ 2w) , the relevant length scale for magnetic field pen- 
etration is the Pearl length A — /dM^ Large current 
fiows near the edges of a wide strip with A/w <C 1^"— 
whereas the distribution of the sheet current for a nar- 



row strip A/w ^ 1 is simply proportional to the distance 
from the center There are, however, no exact analytical 
expressions for strips of arbitrary width in perpendicular 
fields. Numerical calculations with A > can account for 
the sheet current density in some thin geometries under 
perpendicular applied fields.—"— 

The above-mentioned results have been widely used as 
theoretical frameworks for investigating superconductors, 
ranging from studies of fiux penetration in mesoscopic 
superconductors to the macroscopic characterization of 
bulk superconductors. In the present paper, we extend 
the previous results by considering two-dimensional ar- 
rays of superconducting thin strips exposed to perpen- 
dicular fields, and we offer a theoretical framework for a 
new class of superconducting systems, dc magnetic meta- 
materials. 

Metamaterials have recently attracted considerable at- 
tention, because they can provide effective (or macro- 
scopic) electromagnetic properties not found in nature, 
allowing a new set of applications for controlling electro- 
magnetic field i^i^i These effective electromagnetic prop- 
erties come from the internal structure of the metamate- 
rials. Typically, metamaterials are composed of an array 
of structure whose sizes and separation from each other 
are smaller than the wavelength of the involved electro- 
magnetic field. Tailoring such internal structure is the 
key factor in obtaining the desired effective properties. 
A crowning achievement would be electromagnetic field 
cloaking. One can consider, for example, a spherical 
metamaterial with a concentric spherical hole exposed 
to an incident electromagnetic wave; inside the meta- 
material sphere the total electromagnetic field is zero, 
whereas outside the metamaterial sphere the electromag- 
netic field is undisturbed. Effectively, one can find the 
values of the permittivity and permeability tensors that 
ensure the cloaking of electromagnetic fields from trans- 
formation optics^ conformal mapping^ or by consider- 
ing boundary conditions. However, the exact cloaking 
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of a broadband electromagnetic field may involve unphys- 
ical situations, such as a diverging phase- velocity of light 
or violations of causality-energy conditions.— Moreover, 
extreme values of the permittivity and permeability are 
needed and, in most cases, the use of resonant structures 
results in large losses 

The control of a dc magnetic field, either for cloak- 
ing or for other possible applications, represents the zero 
frequency limit of the general electromagnetic field. In 
this dc field, the electric and magnetic fields are decou- 
pled and only the magnetic permeability of the mate- 
rial is relevant.— However, the control of the magnetic 
field requires metamaterials that should be, in general, 
anisotropic and inhomogeneous (i.e., the permeability 
tensor of the material must be a particular function of 
the position), depending on the desired type of control. 

Arrays of thin superconductors such as dc supercon- 
ducting metamaterials^ - can partially provide such char- 
acteristics, since they are intrinsically anisotropic. In- 
deed, maintaining parallel permeability = /io (where 
fiQ is the vacuum permeability), the perpendicular effec- 
tive permeability ^± can be tuned from to /io by chang- 
ing the geometric parameters of the arrayi^ The first 
experimental demonstration of such dc superconducting 
metamaterials providing < 1 was presented in 

Ref. 30, using arrays of thin Pb films. The physics and 
applications of superconducting metamaterials (includ- 
ing dc superconducting metamaterials) were recently re- 
viewed in Ref. Superconductor- metamaterial hybrids 
have also been proposed as antimagnetsi^ 

In the present work we present analytical expressions 
for the magnetic field distribution, as well as for the effec- 
tive permeability of metamaterials consisting of infinite 
regular two-dimensional (2D) arrays of superconducting 
strips in the complete shielding state with A/w ^ 0. 
Such analytical results are useful because they provide 
general trends in the relation between the effective mag- 
netic properties of a metamaterial and the geometry of 
its constituents. Moreover, we numerically calculate how 
these effective properties change when a nonzero A/w is 
considered in the superconducting strips. 

This paper is organized as follows. After an introduc- 
tion to 2D arrays of superconducting strips in Sec.|lTl we 
theoretically investigate two types of 2D arrays of super- 
conducting strips: a rectangular array in Sec. IIIIl and 
a hexagonal array in Sec. IIVI In Sees. IIIIl and IIVI we 
present an analytical investigation of the field distribu- 
tions and effective permeability for the case of complete 
shielding with A/w — >■ 0, and also a numerical investi- 
gation on the effective permeability for the case of finite 
penetration depth with A/w > 0. In Sec. IVlour theoret- 
ical results are discussed and a brief summary is given. 



II. TWO-DIMENSIONAL ARRAYS OF 
SUPERCONDUCTING STRIPS 

We now describe 2D arrays of superconducting strips 
as basic components of magnetic metamaterials. The 
width of the superconducting strips 2w is much larger 
than their thickness d, and their length is infinite along 
the z axis. The wide surfaces of superconducting strips 
are parallel to the xz plane. We consider the thin-strip 
limit case, e — d/2 ^ 0, and hereafter we regard e 
to be an infinitesimal. 

The relation between the local (or microscopic) mag- 
netic field H and the local magnetic induction B is given 
by B = fioH. The macroscopic fields (H) and (B) are 
obtained by respectively averaging H and B in the unit 
cell of the strip array. As shown in Refs. IH and [s^, to 
maintain consistency with the Maxwell equations (H) is 
calculated as the averaged line integral of H, whereas 
(B) is calculated as the averaged surface integral of B in 
the unit cell. Because of the different definitions of aver- 
aging procedure to obtain macroscopic fields, we gener- 
ally have {B) ^ ^o{H), even though B — ^qH holds. 

The macroscopic magnetic response of the 2D arrays 
of superconducting strips is characterized by the rela- 
tion between the macroscopic magnetic field {H) = 
{Hx)x + {Hy)y and the macroscopic magnetic induction 
{B) = {B^)x + {By)y, as (Ba) = ^0/3(^/3), where the 
components of the permeability tensor are fj,xx = 
fiyy = fi±, and fiai3 — for a ^ (3. In other words, 

{Bx) = ii\\{Hx) and {By) = n^{Hy). (1) 

When a 2D array of superconducting strips is exposed 
to a parallel magnetic field (i.e, one along the x axis), 
the field is not disturbed by thin strips of e 0. The 
effective permeability for the parallel field is, therefore, 
equal to the vacuum permeability //|| = /xq for the thin 
strip limit. 

On the other hand, when a 2D array of superconduct- 
ing strips is exposed to a perpendicular magnetic field 
(one along the y axis), the strips strongly disturb the 
field. Because of the magnetic shielding in superconduct- 
ing strips, the effective permeability for the perpendicular 
field is smaller than the vacuum permeability, /ij^ < /io, 
depending on the geometrical parameters for the 2D ar- 
ray, 2a, 26, and 2w. 

Hereafter we study field distributions and effective per- 
pendicular permeability /ij_ for the case where 2D arrays 
of superconducting strips are exposed to perpendicular 
magnetic fields. 

III. RECTANGULAR ARRAYS OF 
SUPERCONDUCTING STRIPS 

Figure [1] shows a schematic of the rectangular array of 
superconducting strips investigated in this section. Su- 
perconducting strips are regularly arranged with a unit 
cell of 2a x 26 in the xy plane; the lattice constant along 
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FIG. 1: (Color online) Cross section of a two-dimensional 
array of superconducting strips in the xy plane. In the nth 
layer at j/ = 2nh, the mth strip is situated at \x — 2ma\ < w, 
where m = 0, ±1, ±2, . . . ± oo and n = 0, ±1, ±2, . . . ± cxd. 



the X axis is 2a, and that along the y axis is 2b, where 
2a > 2w > d and 2b > d. 



A. Complete shielding state, A/w — ^ 

Here we consider field distributions and perpendicular 
permeability of a rectangular array of superconducting 
strips for the complete shielding state (i.e., A/w — ^ 0). 
Such complete shielding can also be seen for a weak-field 
or large- Jc limit (i.e.. Ha/ Jed — )■ 0) in the critical state 
model, where Ha is the applied magnetic field and Jc is 
the critical current density. 



1. Field distribution 

The two-dimensional local magnetic field H = 
Hx{x, y)x + Hy{x, y)y is analytically investigated via the 
complex field 'H(C) = Hy(x, y) -\- iHx{x, y) as an analytic 
function of ^ = a; -|- iy.^^ 

We use the conformal mapping from the C — ^ ~^ W 
plane to rj = u + iv plane^^ 



V = Vi-iO = sn(C/cr,fcr), 



(2) 



where sn(M, k) is the sine amplitude (the Jacobi sn func- 
tion). [See Eq. dMI-] The modulus is the function 
of b/a, and is obtained by solving 



K(fc;.) 



(3) 



a K(fcr)' 

where K(fc) is the complete elliptic integral of the first 
kind^*? [see Eq. (M^] and k'^. = ^1 - k^. The Cr in Eq. (g]) 
is then given by 



Cr = a/K(fc,)=6/K(fc;.). 



(4) 



The upper half of the unit cell in the C plane, |a;| < a 
and < y < b, is mapped onto the upper half plane in 
the 7] plane, < u < ooj^ From Eqs. (jA8p - (|All| . we see 
that the points C = 0, a, a -I- ib, and ib in the C plane 
are mapped onto rj — 0, 1, l/fcr, and oo in the rj plane, 
respectively. 

Using the conformal mapping technique, we obtain the 
complex field for the rectangular array of superconduct- 
ing strips in the complete shielding state, 



no = Ho 



VriC) 



(5) 



Note that the coefficient Hq = Hy{0, b) in Eq. ^ is the 
magnetic field at (x, y) = (0,6). For 6/a -> 0, we have 
fcr -J> 1 and 77r(C) tanh(7rC/2&) from Eq. which 
corresponds to the case of an infinite stack of supercon- 
ducting strips.^ For b/a — )■ oo, on the other hand, we have 
fcr — and r]i-{C) -> sin(7rC/2a) from Eq. (|A7|) . which 
corresponds to the case of an infinite array of coplanar 
superconducting strips.^ 

The perpendicular field component Hy{x,0) = 
Re'H{x) and the sheet current Kz[x) = Hx{x,—e) — 
Hx{x,e) = Im[H(a; — ie) — 'H{x + ie)] in the unit cell 
—a < X < a are obtained from Eq. ([5]): 







for \x\ < 



Hy{x,0) = 



Hori^{x) 



^/r],.{xy - ri,.{wy 
2HaT],{x) 



for w < \x\ < a, 



for |x| < w, 



(6) 



for w < < a. 



where r],-{x) — sn(x/cr, fcr). 

The complex potential defined by ^(C) = jluiOdC,' 
is calculated from Eq. ([5]), and is given by 

CrHa f°° udu 



JvriO y'(u2„l)(y2_fc-2)(y2_^2) 



zF arcsini 



, '^r , (8) 



where F{lp, k) is the elliptic integral of the first kindi^ 
[See Eq. (jAip .] The parameters 7r and in Eq. ^ are 
defined by 



7r = 7]r{w) ^ sn{w/cr,h), (9) 

fci-Cd (w/Cr, fcr), (10) 



1-72 



where cd(u, k) = cn(u, fc)/dn(M, k) is the Jacobi cd func- 
tion (see Appendix |A|). Figure [2] shows the magnetic field 
lines corresponding to the contour lines of Re Q (x+iy) ob- 
tained from Eq. 1^ . The magnetic field concentrates near 
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h/a = 0.5, and (c) h/a — 0.2. Thick horizontal lines correspond to the cross sections of superconducting strips. 



the gaps between the edges of superconducting strips, 
and this field concentration is severe for small stack spac- 
ings 6/a <C 1, as in Fig. [Ifc). 



2. Macroscopic fields 

The local magnetic field H — Hx{x, y)x + Hy{x, y)y is 
obtained from Eq. ([5]), and the local magnetic induction 
is given by J5 = l^-oH. Here we consider the macro- 
scopic magnetic field (H) and the macroscopic magnetic 
induction (B). Because we are interested in the case 
where the apphed magnetic field is parallel to the y axis, 
only y components of the averaged fields (Hy) and the 
macroscopic magnetic induction {By) are relevant. The 
effective permeability is then fj,± — {By) / {Hy) . 

The macroscopic magnetic induction is calculated from 
the averaged surface integral of the local magnetic induc- 
tion in the unit cellj^i^ We therefore define {By) as the 
averaged surface integral of By = fioHy in the region of 
—a < X < a and < z < ai y — b: 



{By 



1 



2aL 
Mo 



dx I dzBy(x,b) 







2a 



z J— a 
+a 

Hy{x, b)dx. 



(11) 



The volume integral of V • S = in t/(?/o), where F(?/o) 
denotes the region of —a < x < a, yo < y < b^ and 
< z < Lz, reduces to 

0=/ VBdV^ [ BdS 

dx dz[By{x,b) ~ By{x,yo)], (12) 

-a Jo 

where the surface integrals of B^ at a; = ±a make no con- 
tribution because Bx{±a,y) — 0. Substituting Eq. (fT2|) 
into Eq. (|lip shows that 



^^y^ " ^ Hy{x,yo)dx 



(13) 



holds for any j/o; that is, the right-hand side of Eq. (fT3|) 
is independent of j/o- 

The macroscopic magnetic field is calculated from the 
averaged line integral of the local magnetic field in the 
unit celli^i^ We therefore define {Hy) as the averaged 
line integral of Hy in the region oi ^b < y < b aX x = a: 



1 

^y) ^2b J ^y'^^'y^'^y- 



(14) 



The surface integral of V x = J in Sz{xo), where 
52(2:0) denotes the region of a;o < 2; < a and —b<y<b 
at z = 0, reduces to 

f Kz{x)dx = I (y xH)-dS 



+b 



-b 



[Hy{a,y) - Hy{xo,y)]dy, (15) 



where the line integrals of H^ at y — ±b make no con- 
tribution because Hx{x,±b) — 0. Equation leads 
to 



1 C^"' f^^ 



Aab 



+a 



dx 



Hy{x,y)dy+ Kz{x')dx' 



Substituting Eq. ^ into Eq. ^ yields 

1 



{H 



1 /■+^ 1 



1 



{B,, 



[By)dy~\ 



/io ' 4a6 j_a 

Equation pT)) reduces to 

{By)/^,0 = {Hy)+{My), 



(16) 



dx / dx'Kzix') 

J X 

{x +a)Kz{x')dx'. (17) 



(18) 
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where we used J^^ Kz{x)dx = 0, and {My) is the mag- 
netization defined by 



{My 



Aab 



xKz{x)dx. 



(19) 



Thus, we have confirmed that the definitions of the 
macroscopic fields given by Eqs. (ITT|) and (|T4)) are con- 
sistent with the relation between the macroscopic fields 
given by Eq. (fT8|) . 



3. Perpendicular permeability 

The macroscopic magnetic induction is calculated by 
substituting Eq. (O into Eq. (fTTj) : 



1 1 /■'^ 1 ra+ib 

— {By) = - / Hy{x,b)dX = - / 

Mo a Jo aj,i, 



nodc 



krO, 



udu 



y^(u2-fcr^)(u2-l)(u2_^2) 



:KK), 



(20) 



where 7r and defined by Eqs. (jH]) and (|10p . respec- 

tively. The macroscopic magnetic field is calculated by 
substituting Eq. (O into Eq. p4)) : 



Hy{a,y)dy 



1 



H(C)dC 



udu 



^(fcr'-u2)(M2_ l)(u2_^2) 



=K«), 



(21) 



where = ^1 - k2 = ^(1 - fc2)/(l - fc272)_ The effec- 
tive permeability /i^j. = {By) / {Hy) is therefore obtained 
as 



^0 aK«)' 



(22) 



When the width of superconducting strips is small 
(ui/cr <C 1), we have n^. ~ k^[l — (1 — k^)w'^ /2c^], and 
Eq. ([22]) reduces to 



Mo 



1 



4a& 



(23) 



When the gaps between the edges of the supercon- 
ducting strips are small {1 — w/a <^ 1), we have Kr — 
krK{kr)il - w/a) and K«)/K(k,) ~ {2/n)\n{i/ k,), 
and Eq. (22) reduces to 



M-Lr 

Mo 



7r5 
2a 



In 



fcrK(fcr)(l-w/a) 



(24) 



0.5 




rectangular array 
h/a = 5,2, 1,0.5. 0.2. andO.l 



0.5 1 
l-w/a 

FIG. 3: (Color online) Effective permeability for perpendicu- 
lar field /ixr as function of 1 — w/a for b/a — 5, 2, 1, 0.5, 0.2, 
and 0.1. The dashed fine is fi±,./fiQ = 1 for b/a oo, and the 
dot-dash line is /ixr//io = 1 — w/a for b/a 0. 



The inverse of /iJ-r/Mo logarithmically diverges when 
1 — w/a — >■ 0. This sharp change in fj.±i- for small gaps 
becomes gradual when we take finite A/w into account, 
as shown in the numerical results in the next subsection. 

When the stack spacings between the wide surfaces of 
the superconducting strips are large {b/a > 2), we have 
kr ~ 4exp(— 7r5/2a) ^ 1, Kr ~ A:r cos(7rw/2a) <C 1, and 
Cr — 2a/7r, and Eq. ([^ reduces to 



M-Lr 

Mo 



1 - 



2a / 
— m cos 

TTO V 



2^7/ 



(25) 



When the stack spacings between the wide surfaces of 
the superconducting strips are small (b/a ^ 1), we have 
fci- ~ 1 — 8 exp(— 7ra/26) and Kj- — 1 — 2 exp[— (7ra/6)(l — 
w/a)], and Eq. (|22|) reduces to 



M-Lr 

Mo 



1 



w 
a 



2b 



In 2, 



(26) 



which is valid for a wide range of w/a. The linear behav- 
ior, /x_Lr/Mo — > 1 — w/a, occurs because the rectangular 
array of strips with b/a ^ corresponds to supercon- 
ducting slabs of width 2w^ Note that Eq. (f26| is invalid 
for w/a <C 1 or 1 — w/a ^ 1. When b/a ^ 1 and 
1 — w/a 1, Eq. (IM)) further reduces to 



M-Lr 

Mo 



irb 
2a 



In 



7r(a - 



w) 



Figure 13] shows a plot of Mir/Mo vs 1 
from Eqs. 0, (Uni), and ([22]). 



(27) 

w/a obtained 



B. Finite penetration depth, A/w > 

We have performed numerical calculations for the re- 
alistic case taking nonzero A/w into account. Our cal- 
culation of the sheet current in superconducting strips is 
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FIG. 4: (Color online) Comparison between the numerical results of the effective permeability /ixr in a rectangular array 
for A/w > (symbols) and the analytical results for A/w — ^ (lines), (a) Numerical results (dots) for A/w — 0.01 and 
analytical results (lines) for A/w — >■ of /^xr//io vs 1 — w/a for b/a = 5, 2, 1, 0.5, and 0.2. (b) Numerical results (symbols) 
for A/w = 0.1, 0.05, 0.01 and analytical results (line) for A/w — >■ of /ixr//io vs 1 — w/a for b/a = 0.2. (c) Numerical results 
(dots) and analytical results (dashed lines) for A/w — >■ of ^xr/A^o vs A/w for b/a = 0.2. 



based on the magnetic energy minimization technique, as 
described in Ref. [H. 

Figure |4] shows a comparison between the numerical 
results of fi±i- for A/w > and the analytical results 
for A/w — > 0. When A/w — >■ 0, the numerical results 
coincide with the analytical ones within the numerical 
accuracy. As shown in Figs.lUJa) andlUJb), /i_Lr increases 
with increasing A/w for any given values of b/a and 1 — 
w/a, similar to what was found in Ref. [2^. When the 
gaps between the edges of the superconducting strips are 
small (1 — w/a < 0.1), we see a pronounced difference 
between ^±r/fJ'0 for A/w > and that for A/w — 0. 
We see a sharp drop in /i^r/^o — >Oasl — w/a->0 
for A/w — )■ [see analytical expressions, Eqs. ([M)) and 
(p7)l ]. in contrast to the slow decrease in //.li/mo with 
decreasing 1 — w/a for A/w — 0.01, because of the finite 
penetration of the magnetic field from the strip edges. 
Figure [4{c) shows /i_Lr as functions of A/w: the for 
A/w < 10~^ agree well with those for A/w — 0, shown 
as dashed lines. We clearly see that for A/w > 0.1 is 
much larger than that for A/w especially for small 
1 — w/a. 
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FIG. 5: (Color online) Cross section of a hexagonal array of 
superconducting strips in the xy plane. In the even layer at 
y = 4nb, the mth strip is at \x — 2ma\ < w, whereas in the odd 
layer a,t y — (4n-|-2)6, the mth strip is at \x~ {2m + l)a\ < w, 
where m = 0, ±1, ±2, . . . ± oo and n = 0, ±1, ±2 . . . ± oo. 



1. Field distribution 



IV. HEXAGONAL ARRAY OF 
SUPERCONDUCTING STRIPS 

In this section, we consider a hexagonal array of su- 
perconducting strips shown in Fig. [51 and compare its re- 
sponse with that of the rectangular array shown in Fig. [1] 



A. Complete shielding state, A/w 

In this subsection we consider field distributions and 
perpendicular permeability of a hexagonal array of su- 
perconducting strips for A/w 0. 



Here we use the conformal mapping from the C = x+iy 
plane to the rj — u + iv plane?^ 



V = VhiO = sn(C/ch,fch)- 
The modulus fch is obtained by solving 



2b 

a 



K(fch)' 



(28) 



(29) 



where = ^Jl — k"^. The relation between fch from 
Eq. (gni) and k, from Eq. dS]) is given by fch = (1-A:;.)/(1 + 

K) = (l-V^l - fc?)/(l + v/l - ^?)- ch in Eq. dMl) is given 
by 



Ch = a/K(fch) = 2b/K{k'^). 



(30) 
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FIG. 6; (Color online) Magnetic field lines in a hexagonal array of superconducting strips for w/a — 0.8: (a) b/a — 1, (b) 
b/a — 0.5, and (c) b/a = 0.2. Thick horizontal lines correspond to the cross sections of superconducting strips. 



The upper half of the unit cell in the C plane, |a;| < a and 
< y < 26, is mapped onto the upper half plane in the 
T] plane, < v < ooM- The points C = 0, a, a + 2ib, and 
2ib in the ^ plane are mapped onto ?/ = 0, 1, 1/fch, and 
oo in the rj plane, respectively. 

The complex field H(C) = Hy{x, y) + iHx{x, y) describ- 
ing the local field distribution for a hexagonal array in 
the complete shielding state is given by 



m) = H, 



where 7h 



r]u{w) = sn(w/ch,fch), and /3h 



(31) 



r]h{a 



w + 2tb) = ^{k-^^ - 7^)7(1 - il) = (l/fch)dc(«;/ch, fch). 
(See Appendix A for Jacobi dc function.) Note that 
the coefficient Hq = Hy{0,2b) in Eq. (|3T1) is the mag- 
netic field at {x,y) = (0,26). For b/a — )■ 0, we have 
fch -> 1 and 77h(C) ~^ tanh(7rC/46), which corresponds to 
the case of an infinite stack of superconducting stripsJ' 
For 6/a — > oo, wc have fch — > and rjh{C) — >■ sin(7r^/2a), 
which corresponds to the case of an infinite array of copla- 
nar superconducting strips.^ 

The perpendicular field component Hy{x,0) = 
KeH{x) and the sheet current Kz{x) = H^{x,—e) — 
Hx{x, e) = Im \H{x — ie) — H^x + ie)] are obtained from 
Eq. ([5T|). and are given by Eqs. ^ and ([7]) by replacing 
rir{x) by 7?h(a;) = sn(w/ch, fch)- 

The complex potential defined by G{C) = /2^if, ^(C')'^C' 
is calculated from Eq. 



udu 



fch 7^„(C) ^{u^ - 1)(m2 



ChHo 



zF arcsin < 



9 J 

7h 



1(32) 



where k'^ = -^/l — and /th is defined by 
Kh — -1 / r — n 7^ ?:7^ — 1 



fch -l + (l-7if) 



2^2 



fc: 



-1-1/2 



cn"* {w/ch, fch) 



(33) 



Figure ini shows the magnetic field lines that correspond to 
the contour lines of ReQ{x + iy) obtained from Eq. (|32l) . 
When the stack spacing is large {b/a > 1), field distri- 
butions in a hexagonal array are similar to those in a 
rectangular array, as seen in Figs, ^a) andlHl^a). When 
the stack spacing is small {b/a < 1), on the other hand, 
we see a striking difference between the rectangular and 
hexagonal arrays. The magnetic field concentrates only 
near the gaps between the edges of the strips in a rectan- 
gular array [Fig.[21[c)], whereas the magnetic field is large 
in most of the region in a hexagonal array [Fig. |6ljc)]. 



2. Macroscopic fields 

The definitions of the macroscopic magnetic induction 
(By) and the magnetization (My) for a hexagonal array 
arc the same as those of a rectangular array, and are 
given by Eqs. ([TT|) and ([TO)) , respectively. 

The definition of the macroscopic magnetic field (Hy) 
given by Eq. (jl4p for a rectangular array, on the other 
hand, must be modified so that (Hy) for the hexagonal 
array is consistent with the macroscopic relation given 
by Eq. ^TE\\ . We therefore define (Hy) for a hexagonal 
array as 



I Hy{a,y)dy- / Hx{x,2b - e)dx 
Jo Jo 



(34) 



3. Perpendicular permeability 



The resulting effective permeability for hexagonal ar- 
ray /z_Lh = {By)/{Hy) is given by 



/^-Lh _ 26 K(Kh) 
Mo a K{k'J ' 

where Kh is defined by Eq. 



(35) 
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FIG. 7: (Color online) Effective permeability for perpendicu- 
lar field iJ.±h as function of 1 — w/a for b/a = 5, 2, 1, 0.5, 0.2, 
and 0.1. The dashed line is /xxh/A^o = 1 for b/a — >■ oo. The 
dot-dash line of /ix//io ~ 1 — m/a is shown for comparison 
with Fig. 



When the width of superconducting strips is small 
{w/ch ^ 1), we have Kh — fch[l — (1 — fch)^^/2ch], and 
Eq. §51 reduces to Eq. (|23l) . 

When the gaps between the edges of the superconduct- 
ing strips are small (1 — w/a ^ 1), we have 



Mo 



2a 



In 



[fcf^(l-fc2)]i/4K(fcO(l-W«) 



Equation is further simplified for b/a<^ 1 as 



Mo 



7r6 
2a 



7ra 

Trr +ln 



46 



7r(a — w) 



(36) 



(37) 



When the stack spacings between the wide surfaces of 
the superconducting strips are large {b/a > 2), we have 
fch — 4exp(— 7r6/a) ^ 1, Kh ~ fch cos^(7r?i;/2a) <C 1, and 
Ch — 2a /it, and Eq. ((35|) reduces to Eq. ((25|) . In other 
words, the magnetic response of the hexagonal array is 
almost the same as that of the rectangular array for large 
b/a, as expected. 

Figure [7] shows a plot of /i±h/Mo vs 1 — w/a obtained 
from Eqs. and Note that for a given b/a, 

/i_Lh of the hexagonal array shown in Fig. [5] is smaller 
than /xxh of the rectangular array shown in Fig. [1] espe- 
cially when b/a ^ 1. The effective permeability of the 
hexagonal array is very small and slightly dependent on 
the gaps between the edges of the superconducting strips; 
M±1i/Mo ^ 1 when the stack spacing is small {b/a <^ 1) 
for the wide range of w/a, < 1 — w/a ^ 0.5. The hexag- 
onal array is, therefore, more advantageous in obtaining 
highly anisotropic permeability than the rectangular ar- 
ray. 



B. Finite penetration depth, A/to > 

Figure [5] shows a comparison between the numerical 
results of /i_Lh for A/w > and the analytical results for 
A/w — ^ 0. The overall tendency in the difference between 
fi±h for A/w > and that for A/w — ^ shown in Fig. [S] 
is similar to that in Fig. 21 It is interesting to compare 
fi±h of a hexagonal array and /i_Lr of a rectangular array 
for small stack spacings {b/a <^ 1), in that the /ij^h for 
A/w > is close to that for A/w — ?► 0, even for small 
1 — w/a. 



V. DISCUSSION AND SUMMARY 

In this paper, we have investigated anisotropic effec- 
tive permeability of infinite arrays of thin superconduct- 
ing strips, which can be regarded as essential structures 
of a dc magnetic metamaterial. When an applied mag- 
netic field is parallel to the superconducting strips, the 
effective permeability is given by /iy = /iq. When an 
applied magnetic field is perpendicular to the strip, on 
the other hand, the effective permeability is ii± < fiQ, 
because of the magnetic shielding in the superconduct- 
ing strips. The perpendicular permeability fi± becomes 
small when the stack spacings are small {b/a <C 1), or 
when the gaps between the edges of the strips are small 
{1 — w/a ^ 1). We have investigated two types of 2D 
arrays: the rectangular array shown in Fig. [1] and the 
hexagonal array shown in Fig. [S] The hexagonal array 
is more advantageous for obtaining large anisotropy of 
the effective permeability than the rectangular array, be- 
cause /ij^h is much smaller than /ij^i- when 6/a ^ 1 and 
l-w/a<0.5. 

To realize magnetic cloaking or other possible applica- 
tions for controlling magnetic fields, magnetic metama- 
terials with large anisotropy (/i±/Mo ^ 1) are needed. 
As seen in Figs. [3] and [3 m±/Mo ^ 1 can be achieved by 
means of an array of superconducting strips with small 
gaps 1 — w/a ^ 1 and small stack spacing b/a ^ 1. For 
finite A/w, however, /i^//io ^ 1 is difficult to realize, be- 
cause of field penetration near the edges of the supercon- 
ducting strips. Wide superconducting strips are favor- 
able for obtaining small A/w and small m^/MOi but con- 
trol of magnetic fields for magnetic cloaking, for example, 
requires narrow superconducting strips. We therefore 
need to consider optimization of the strip size and con- 
figuration of the array of superconducting strips. Such 
optimization depends on the details of the application for 
magnetic metamaterials. 
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FIG. 8: (Color online) Comparison between the numerical results of the eflFective permeability fi±h in a hexagonal array for 
A/w > (symbols) and the analytical results for \/w — >■ (lines), (a) Numerical results (dots) for A/w — 0.01 and analytical 
results (lines) for A/w — >■ of i-L±h/fJ'0 vs 1 — w/a for b/a = 5, 2, 1, 0.5, 0.2, and 0.1. (b) Numerical results (symbols) for 
A/w — 0.1, 0.05, 0.01 and analytical results (line) for A/w — >■ of /ixh/Mo vs 1 — w/a for b/a = 0.2. (c) Numerical results 
(dots) and analytical results (dashed lines) for A/w — >■ of ^t±h//io vs A/w for b/a — 0.2. 



Appendix A: Elliptic integrals and functions 



(i.e., delta amplitude) are defined via u — F{ip, k) by 



In this paper, some elliptic integrals and functions ap- 
pear, but those notations and definitions are not unified 
in the literature. In this appendix, we therefore summa- 
rize those expressions to avoid ambiguity. The expres- 
sions below are from the textbook by Gradsytein and 
Ryshik.36 

The elliptic integral of the first kind F{(p, k) and the 
complete elliptic integral of the first kind K(fc) are given 
by 



Fi^,k) 



sm 



dt 



K(fc) = F{-K/2,k) 



^{\-f^){l-kH-^y 
dt 



^{\-t^){i-kH^y 



(Al) 
(A2) 



respectively. 

The Jacobian elliptic functions sn(u, fc) (i.e., sine am- 
plitude), cn(u, fc) (i.e., cosine amplitude), and dn(M, fc) 



sn(u, /c) = sin(/5, 
CT\(u,k) — cos(/3. 



dn(u, k) = y 1 — k^ sin ip 



(A3) 
(A4) 

(A5) 



respectively. We also use cd(u, fc) — cn(u, fc)/dn(u, fc) 
and dc(M, k) = dn(u, fc)/cn(u, k) = l/cd(u, k). 

Here are simple expressions of sn(M, k) for specific k or 

u: 

sn(u, 1) ~ tanh(u), (A-6) 
sn(u, 0) = sin(w). (A7) 
sn(0,/c) = 0, (A8) 
sn(K,fc) = 1, (A9) 
sn(K + iK',/fc) l/k, (AlO) 
sn{iK\k) ^ oo, (All) 



where K = K{k) and K' = K(/c') = K(Vl - fc^). 
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